Non-equilibrium Entanglement and Noise in Coupled Qubits 
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We study charge entanglement in two Coulomb-coupled double quantum dots in thermal equi- 
librium and under stationary non-equilibrium transport conditions. In the transport regime, the 
entanglement exhibits a clear switching threshold and various limits due to suppression of tunneling 
by Quantum Zeno localisation or by an interaction induced energy gap. We also calculate quantum 
noise spectra and discuss the inter-dot current correlation as an indicator of the entanglement in 
transport experiments. 
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Precise engineering and preparation of entangled states 
forms the backbone of many quantum information 
schemes Q. The complete control of interactions be- 
tween two or more parties is a sangraal that is not with- 
out cost. For example, in superconducting nano-circuits 
there has been much success in devising schemes for 
tunable capacitative coupHngs but thermal fluctua- 
tions, background noise and limited control over natural 
interactions must be dealt with and overcome in increas- 
ingly imaginative ways. 

In this Letter, we take a slightly different point of 
view and ask for the degree of entanglement between 
two parallel, interacting electronic conductors under the 
'un-favourable' condition of stationary currents passing 
through both of them. As this is clearly a mixed-state 
situation, we specifically consider a non-equilibrium ver- 
sion of the concurrence as entanglement measure for an 
electron charge double qubit (DO), realised in Coulomb- 
coupled double quantum dots [^[g that are strongly cou- 
pled to external electron reservoirs at high voltage bias. 
We compare this to the same closed device in equilib- 
rium with a heat bath, and our findings suggest that 
such an approach, although probably not directly rel- 
evant for quantum information purposes, sheds a new 
light on the relation between entanglement and the elec- 
tronic transport process itself. In particular, effects like 
suppression of non-resonant tunneling and the quantum 
Zeno effect (QZE) have a direct impact on the entangle- 
ment, to which we also establish a further link by cal- 
culating the non-equilibrium quantum shot-noise tensor 
whose off-diagonal elements, as a function of the system 
parameters, show a behavior very similar to the concur- 
rence. 

Model.- For the sake of clarity, we define the double 
qubit by 'left' and 'right' orbital charge states \ai), a 
L,R oi one additional electron on top of the many-body 
ground state |0i) (limit of intradot Coulomb interaction 
Uin — > oo) of two double quantum dots i = 1,2 which 
are coupled by a single matrix element U for inter-Aoi 




FIG. 1: Transport double qubit model: left-right charge states 
in two Coulomb-coupled double quantum dots with one ad- 
ditional electron each and 'on-site' {LL, RR) interaction U, 
coherent tunnel couplings Ti and T2 and electron reservoir 



tunnel rates F 



1/2 

L/R- 



'same site' interactions (left- left and right-right), cf. Fig. 
([T|l. Tunnelling of electrons occurs only within but not 
between the qubits due to coupling Tj in each double 
dot. Using projectors onto these orbital states, (n^''' = 
^LH = \L){R\i,...), the total Hamiltonian is 



\L){LU, n) 
Ho 



E 
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The electron spin label is suppressed here and in the fol- 
lowing, as only charge states (acting as pseudo-spin) play 
a role. This description has turned out to be useful for 
modelling charge-related properties such as decoherence 
and noise in individual double quantum dots 0, 0] . 

We 'open' the DQ by coupling it to four external elec- 
tron reservoirs, H ~ Hq + Ht + 'H^cs, with Tires = 



1,2 Z^aeL,R 



{a = L/R refers to left 



1,2) and 
with Hub- 



and right reservoirs for qubit number i, i = 

nr = E,=i,2EoeL,flE.(n"'cL^L + ^-c.), 

bard operators = |Oi)(Q!i| that couple the qubits to 
the continuum. 
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FIG. 2: Left: Grey-scale plot of double qubit equilibrium 
concurrence as a function of interaction U /Tc and tempera- 
ture T/Tc (white is zero, black is maximum C = 1 ). In all 
of the following results both devices have identical parame- 
ters r*^'^ = Fl, = Tr, Ti = n, for i = 1,2. Right: 
Concurrence of non-equilibrium double qubit as a function 
of interaction U /Tc and reservoir tunnel rate r_R for large 
Vl/Tc — 50. Zero entanglement occurs below a threshold 
oc I/Tr in the weak tunnelling regime, and for very strong 
tunnelling Tr ^ Tc due to Zeno-localisation, cf. text. (Black 
is the maximum C = 0.3). 



Equilibrium entanglement. - If the DQ is disconnected 
from the reservoirs (TYt = 0) but in contact with a heat 
bath at temperature T — 1/ [3^ the equilibrium entan- 
glement between qubit 1 and 2 is easily obtained from 
the concurrence (a well-known entanglement measure 
of mixed states of two qubits) C(/3) for the canonical 



ensemble state p(/3) 



/Z, Z = Tre-'^^o. The 



eigenvectors of Tin correspond to eigenvalues i?o = 0, 
El = [/, and E± = {U ± ^IQT^ + m)/2 and are ex- 
pressed in the basis of singlet and triplet states, 5*0 = 
1/V2(|iii?2> - |i?ii2», T+ = IL1L2), T_ = \R1R2), 
To = l/V2{\LiR2) + \R1L2)). For simplicity we re- 
strict ourselves to the unbiased, symmetric case Ei = 0, 
Ti = Tc(i = 1,2). It turns out that the equilibrium case 
already exhibits some interesting features, cf. Fig.©. At 
any finite temperature T, the entanglement is zero below 
a certain threshold value of the interaction U where the 
state is too mixed in order to be entangled which is, 
e.g., in analogy with the corresponding transition in the 
(abstract) example of the Werner state 0|. Furthermore, 
the concurrence shows a non-monotonic behaviour as a 
function of U at fixed T, with an entanglement maximum 
at an optimal J7-value. 

Stationary transport.- The limit i ^ 00 in the dy- 
namical evolution of the reduced DQ density operator 
p defines a stationary non-equilibrium state p^o which 
usually is much more difficult to determine than in the 
equilibrium case. Transport properties of models like 
Eq. can be analysed by using various non-equilibrium 
techniques. Here, we consider a specific limit of infinite 
source-drain bias in order to obtain quasi-analytic results 
from a generalised Master equation, p = L[p]. The super- 



operator L is parametrized by the Markovian DQ-lead 
tunnel rates F^ = 27r^^ l^fc*1<5(£ — ^kia) (of which the 
energy-dependence is neglected), and the DQ parame- 
ters Ei = 0, Ti = Tc{i — 1,2). Analytical expressions 
for the stationary solution of the 25 coupled equations of 
motion (EOM) can then be found in an approximation 
where the broadening due to tunnelling of the DQ levels 
is neglected, which for Si = 0, however, is only a very 
crude approximation. 

One obtains better results for the stationary currents 
{Ii)oc by second order perturbation theory in the intra- 
dot tunnel couplings Ti, which clearly show a tunnel- 
broadened resonance 

(— e is the electron charge). In this limit, the resonance 
is determined by the energy gap U between the localised 
eigenstates of the DQ: at large U, the triplet T+ = \L1L2) 
becomes populated (note that this state is always avail- 
able because of the infinite voltage approximation). The 
energy gap to any other state involving delocalized elec- 
trons (e.g, the triplet Tq or the singlet Sq) then suppresses 
the elastic current. In analogy to single charge qubits, 
where the energy gap is given by the internal bias e, we 
expect this suppression to be lifted in the presence of 
inelastic processes [lo| . 

Furthermore, as a function of the coupling F^ to 
the drain, the current first increases and then becomes 
smaller again. With the drains acting as broadband mea- 
suring devices (electron on right side or not), strong cou- 
plings F^ — > 00 completely freeze the charges on the left 
sides which is a 'transport version' example of the 
QZE. Alternatively, this localisation can be interpreted as 
an infinite level broadening and the corresponding sup- 
pression of the local spectral density due to the decay 
to the drain. Finally, the behaviour of the current, cf. 
Eq.(|2J), follows the occupation of the entangled singlet 
state as illustrated in the inset of Fig.JJ)). The main 
current contribution therefore stems from two-particle 
tunneling events, which in turn motivates our later com- 
parison of the concurrence with the current fluctuations. 

Non-equilibrium entanglement.- We now define the 
non-equilibrium entanglement via the concurrence C of 
the stationary state Ppoc, where P is the projection onto 
doubly occupied states including proper normalisation; 
i.e. we calculate the concurrence when both double dots 
have a single electron in them and there are thus two 
two-state systems to be entangled. The projection P cor- 
responds to taking the limit F^ — > 00 where both qubits 
are always occupied with one single electron. For exam- 
ple, for [/ = and F^ ^ cxd, the stationary state of a 
single charge qubit is described by the (Bloch) vector of 
pseudo-spin Pauli matrices {e = Sl ~ £b) 



( 2T,e TrT, Fy4 + 



(3) 
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FIG. 3: Main: Diagonal noise spectrum F{uj)i^i — 
S{u)\^\/2el . The resonance at = 2Tc splits into new res- 
onance points at the Bohr Frequencies A± = Ei — E± — 
1/2(C/ =F ^/T6Tf~nP). Top right inset: Stationary cur- 
rent Jstat- Top left inset: Occupation of several singlet 
and triplet states. Bottom Inset: The cross-correlation fre- 
quency spectrum F{lj)i^2 = S{uj)i^2/2el . Again resonance 
points manifest, however the correlation is always zero for 
U = and can assume negative values for U 0. 



with TV = + + 2T^ and in the L-R basis where 

az = \L){L\ — \R){R\ etc. For U 7^ 0, we numerically 
checked that Ppoc = linir^^cxi Poo which means that 
the two-qubit concurrence C defined in this way does 
no longer depend on the left tunnel rates. This is a good 
description of non-equilibrium entanglement in a real sys- 
tem as long as T]^ ^max([/, F^, Tc, e). 

As C is zero to second order in Tc, we use numerical 
results, cf. Fig.(|2J), which shows an intriguing behaviour 
of the concurrence as a function of U and the tunnel 
rate = Ffl. We find a switching threshold in that 
below an interaction strength U ~ 2T'^ /Tn the entan- 
glement is zero: for small Tr, the stationary currents 
become very small, cf. Eq. (j^J, and thus strong interac- 
tions are required in order to entangle the dots. The DQ 
state becomes strongly mixed for Tr (note that we 
have not taken into account any additional, internal re- 
laxation processes here) ; its zero entanglement along the 
axis F/j = is in fact a continuation of the point U = 
where the states of both qubits are located at the origins 
of their Bloch spheres, cf. Eq. (j^J. 

On the other hand, for very large F/j one runs again 
into the QZE with electrons becoming trapped on the left 
{{az) — *■ 1, cf. Eq. 0), and Ppoo approaching the (pure) 
localised state \L1L2) which has zero entanglement. Fi- 
nally, an increase from small to larger F/j at fixed U yields 
the re-entrance behaviour visible in the 'teardrop'-shaped 
region of large entanglement in Fig. 10) . 

Non-equilibrium noise: formalism.- Turning now to 
our description of non-equilibrium shot-noise and its rela- 
tion to entanglement, the stationary state poo on its own 



is not sufficient in order to describe intrinsic properties 
of the DQ: for example, only limited information on the 
spectrum can be obtained from stationary quantities like 
the current. In contrast, the shot-noise spectrum exhibits 
resonances at the transition frequencies of the system and 
contains furthermore useful information on its relaxation 
and dephasing properties 0, 0, 0, We will now 

also show an emergent resemblance in the behaviour of 
the current cross noise and the non-equilibrium concur- 
rence as a function of the system parameters. 

In general, the finite-frequency noise has contributions 
from particle currents as well as contributions from dis- 
placement currents In our case {T\^ :$> F^), how- 
ever, it is a good approximation to consider only particle 
currents. Our starting point is the generating function 



G(si,. 



E 



(ni),...,(nm) 



(t) (4) 



which, for an arbitrary number of m qubits, contains 
the complete information on the tunnelling process as 
a function of time via the counting variables s = {si} 
and the conditional density matrices - (^^■j fQj. 

Ui tunnelling events ('jumps') to the drain i after time 
t. In matrix form, the EOM of the generating function 
follows from the Liouville equation for the conditional 
density matrices and reads G(s,t) ~ M{s)G{s,t) with 
formal solution G(s,t) = exp[iM(s)]G(s, 0). General 
expectation values can be extracted from derivatives of 
Tr[G(s,r)] with respect to the counting variables. In 
particular, the symmetrized noise correlation function 
S{uj)ij = J^^e"^m6h{t + T),SIj{t)}) between qubit 
i and 7 can then be written MacDonald formula 

El El 



S{'^)i,] 



/•OO Q 

/ (irsin(wr) — 
Jo dr 



(5) 



i)ijTr[G(s, r)]|s=i with the differential 



where {niUj) 
operator Dij = dsi^s^ -r 
sion following Flindt et al. by introducing jump op- 
erators Li for qubit sources i and writing dr{ninj 



Sijd^ We simplify this expres- 



(ni),...,(n„) 



(t)] + {i <-> j). This can 



be further evaluated by Laplace transforming the EOM 
dtG{s,t) = {Lq + s.iLi)G'(s, i) and taking deriva- 
tives in counting variables, giving dsiG{s, —i(jj)\s=i = 
Fi^LiF^po, where = {—iuj — L)^^ and po is the 
steady state initial condition. Using the projections 
= -P/iu} - Ruj, Ru> = Q{tuJ + L)-^Q, {P = po<E> 1, 
Q = I — P) with Ppo ~ po and Qpo = leads to 



S{uj)i. 
-2e2 



RcTr 



LjPo 



+ {i^ 



j)-(6) 



In the zero frequency limit, we verify that the noise 
is determined as usual 0, 0| by the lowest eigen- 
value Ao(s) of the matrix M{s), namely by the long- 
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FIG. 4: Left: The cross-correlation zero- frequency spectrum 
Fano Factor -F(0)i,2 = 5(0)i,2/2e7. The resemblence to the 
concurrence, Fig. is qualitative. (White is minimum, 

Fi,2 = —0.12, black is maximum, Fi,2 = 0.84). Right: 
The switching phenomenon in the concurrence is more clearly 
seen, as well as a negative to positive re-emergence in the first 
derivative of the noise around the same point. 



time behaviour G{s,t — * oo) oc exp[iAo(s)] and therefore 
^(0),,, = 2e2A:,Ao(s = 1). 

Non-equilibrium noise: results. - The currents through 
our two parallel charge qubits give rise to two diagonal 
and one ofF-diagonal component in the tensor S{uj)ij of 
the noise spectrum. In Fig. , we present results for the 
diagonal noise, i.e. the noise spectrum S{uj)i^i ~ S{lo)2.2 
of the individual, interacting qubits. This spectrum 
clearly displays resonances at the Bohr frequencies as 
given by the excitation energies of the closed system. At 
[/ = 0, there is one single resonance at = 2Tc that splits 
up when U is increased. Similar to light emission spectra 
in real molecules, frequency-dependent shot-noise spec- 
tra thus provide direct information about the correlated 
energy levels in artificial molecules. 

The cross-noise spectrum exhibits a somewhat more 
complicated resonance structure (inset of Fig.Q). More 
interesting is however the behaviour of the cross- 
correlation Fano factor at zero frequency, F(0)i.2 = 
S'(0)i_2/2e/, which becomes positive as U increases. This 
positive cross-correlation is an indication of correlated 
emission of electron pairs into different exit right leads 
|20j . As a function of [/, cf. in Fig.Q, there is a 
strong analogy between the cross-noise F(0)i.2, its first 
derivative -F"(0)i^2, and the non-equilibrium concurrence 
C, at least on a qualitative level. In particular, the 
non-analytic switching of C with increasing U from un- 
entangled to entangled states translates into a strongly 
delayed (though smooth) onset of the increase in F(0)i.2, 
and the transition of F'(0)i^2 from negative to positive. 
For example, in the figure we see F'(0)i.2 for Tc = 0.5 
becomes positive around [/ = 1 in agreement with the 
switching of C at [/ - 2T'^/Tr = 1. This analogy be- 
tween noise and entanglement so far holds on a qualita- 
tive level only. 

Finally, we mention that we have not included the ef- 



fect of dissipation in our calculations so far. Weak deco- 
herence processes can in principle be easily incorporated 
through additional terms within the master equation. In 
it was shown how to use the resulting changes in the 
noise spectrum in order to extract, e.g., relaxation and 
decoherence times Ti and T2. This can also be done for 
the interacting qubits discussed here. 

Conclusions We have shown how the entanglement of 
a non-equilibrium double qubit differs from its thermal- 
equilibrium relative by exhibiting a l/F switching thresh- 
old for weak tunnelling rates F. The cross-correlation 
noise reflects this threshold at w = and shows reso- 
nances at the Bohr frequencies of the double qubit for 
finite uj. Future theoretical work may include clarifying 
this relationship and checking the influence of decoher- 
ence on the correlated noise power spectrum. 
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